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We review the status of Standard Model predictions for lifetimes and 
mixing rates of charmed mesons. It is shown that the short distance 
approach is able to reproduce t(D + )/t(D°) at leading order in the l/m c 
expansion. SU(3) violating effects from interactions with the soft hadronic 
background are identified as the dominant contribution to the D-D mixing 
rate. We discuss the contribution from operators of dimension nine, which 
is able to enhance the neutral charm width splitting by a factor of order 
ten and comment on possible CP violation in mixing. 
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1 Introduction 



Mixing and CP violation of charmed mesons is being probed with unprecedented 
precision at LHCb pQ. Setting bounds of some 10 3 TeV on the effective scale of various 
AC = 2 operators [2], it severely constrains possible extensions of the Standard 
Model. The D sector is complementary to B and K in offering a handle to probe 
flavour changing neutral currents among weak isospin up quarks [3]. The recent 
evidence for CP violation in D° — > K~K + , 7r~7r + decays [4 J triggered additional 
interest in charm phenomenology. In the Standard Model, charm physics is dominated 
by the first two generations, and CP violation is small. In absence of a Standard 
Model background to interfere with, CP violation in charm is commonly considered 
as a promising possibility for the search of physics beyond. Large penguin effects, 
however, partially compromise this feature. 

Flavour oscillations of neutral mesons arise due to non-zero mass and width split- 
tings, AM and Ar, between the stationary eigenstates. Charm mixing is experimen- 
tally established the most recent HFAG averages [6] for the mixing parameters 
are 

^=(0.63^)%, 

AT (1) 
V= — = (0.75 ±0.12)%. 

Two frameworks are being explored to describe these quantities theoretically: the 
short distance, or inclusive, approach assumes quark-hadron duality and is based 
on parton-level perturbation theory. In the long distance, or exclusive, approach 
the width difference is expressed as a sum over final states common to D° and D 
decays [TJ. In both frameworks, substantial hadronic uncertainties impede definite 
predictions of the Standard Model expectation. 

Our work is anchored in the inclusive framework. A straightforward application 
to charm mixing is known to run into severe trouble, what is usually credited to a 
failure of quark-level perturbation theory at the charm threshold. In the second sec- 
tion below we discuss several issues which we think make it worth to have a closer 
look, though. We present first results of a Standard Model calculation of the D + -D° 
lifetime difference, which avoids the SU(3) flavour interference challenging the calcu- 
lation of mixing rates. We find an intriguing agreement with experiment, somewhat 
clouded by large parametric uncertainties due to missing lattice input. As regards 
mixing, SU(3) breaking from hadron state interactions has been conjectured to be 
the dominant contributions and possibly explain the observed mixing rates within 
the short- distance picture. We report on a calculation of a specific class of such con- 
tributions, and find that we can enhance the D D decay width difference by a factor 
of order 10. As already previously, we argue that from the current theoretical status, 
we can not exclude CP violation in mixing of up to one per cent. 
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2 The D + - D° lifetime difference 

The heavy-quark expansion [8J allows to asses decays of mesons containing one heavy 
quark. It has proved to be an excellent tool to describe mixing and decays in the B 
meson sector. The power of heavy-quark techniques became quite impressive earlier 
this year, when LHCb, CDF, and DO [9] reported on new measurements of the inverse 
lifetime 1/r (B s ) using angular analysis [TD] in B s — > ifxp. Contrasting the most recent 
Standard Model determination of B meson lifetime ratios [TT] with the numbers 
presented from LHCb [9] in Moriond reveals an amazing agreement: 

r(B s ) t(B s ) 

——^ = 1.001 ±0.014 , —^4 = 0.996... 1.000 . (2) 

T ( B d) exp T ( B d)sM 

In the same analysis, LHCb also reported on the first measurement of a non-zero 
width splitting Ar in the B s meson system exceeding 5a [9]: 

AT{B S ) = (0.116 ± 0.019) ps" 1 . (3) 

The current HFAG average [6] is (this average does not yet include the above number, 
but uses the one previously published [12]): 

AT(B S ) = (0.100 ± 0.013) ps" 1 . (4) 

These numbers had been of particular interest, since the B s width difference is be- 
lieved to be most sensitive to violations of quark hadron duality and receives sub- 
stantial contributions from hadronic scale dynamics (lattice bag parameters departing 
from one), perturbative QCD, and subleading 1/rrit, corrections [13J: 

Ar [B s ) = AT (B s ) x (1 + <5 lattice + <5 QCD + <5 HQE ) 

= 0.142 ps" 1 (1-0.14 -0.06 -0.19) . (5) 

Comparing to the most up-to-date Standard Model calculation [TT] . 



f {B J ) r = 0100 ± 0013 = 1.15± 0.32, (6) 
Ar(5 s ) SM 0.087 ±0.021 ' v ; 

shows that, even under more adverse conditions, the heavy-quark expansion managed 
to provide a solid prediction up to (at least) 30% accuracy. 

In charm, our image of the inclusive framework is still rather hazy. We do not 
know whether the radius of convergence in heavy-quark expansion is large enough 
to allow for a perturbative treatment of mass-suppressed corrections down at the 
charm threshold. On top of that, the lower scale brings about substantially richer 
QCD dynamics. Applied to mixing, the naive leading-order predictions in the neutral 
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WA PI 

Figure 1: Weak annihilation (WA) and Pauli interference (PI) diagrams at leading 
order in QCD. They appear in the l/m c expansion at dimension six and are the 
dominant contributions to T (D°) and T {D + ). Shaded circles indicate insertions of a 
AC = 1 operator. 

charm system fall substantially short of the observed values (DQ) for mass and width 
differences. There are some observations, on the other hand, suggesting that it might 
be worth to have a closer look at the inclusive short distance picture. To begin with, 
have a look at what has actually changed compared to the B meson system: heavy- 
quark expansion is a series expansion in hadronic scale over energy released in decay 
modes generating the D-D or B-B transition. The dominant contributions to Ar 
are from the D S D S final state in B s and from KK and ttti in D, respectively. It is 
striking, that the energy releases in D are not significantly smaller than those in B s , 
where the short- distance toolkit works really well: 

B° s D+Dj 1.4 GeV 

D° 7T7T 1.6 GeV 

D° -> KK 0.9 GeV 

With hindsight we also know that the expansion parameter in AT(B S ) turned out to 
be around 1/5, implying that the relevant hadronic scale is significantly below the 
1 GeV it is commonly expected to be. An earlier calculation of subleading corrections 
in charm mixing [H] likewise did not show signs of a breakdown of the perturbative 
approach: the charm width difference receives corrections from next-to-leading order 
QCD at a level of below 50%, and l/m c -corrections of 30%. 

Throughout the last two years, we have advocated [14] that lifetime measurements 
of charmed hadrons might offer a useful way to assess the applicability of heavy-quark 
methods to the charm sector on a more quantitative basis: while the leading orders 
in l/m c describe spectator model quark decays and are (almost) not sensitive to the 
meson's light degrees of freedom, D° mixing and the D + - D° lifetime difference arise 
due to weak interaction with the light valence quarks. They contribute through terms 
of dimension six onwards and test the heavy-quark expansion at the same order in 
l/m c . The D + - D° lifetime difference, generated by weak annihilation and Pauli 
interference (see Fig. [1]), goes back to effects very similar to the ones accounting for 
the neutral charm width splitting. As opposed to mixing rates and CP violation, 
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parameter input rel. error 

Td 206.7 ± 8.9 MeV [T7] 5.48% 

B x 1 ± 1/N C 20.63% 

B 2 1 ± l/N c 0.60% 

ei ± 1/10 49.47% 

e 2 0±1/10 8.81% 

m lGeV...2m c 

/i lGeV...2m c t\lf 8 % 

Aqcd 222 ±27 MeV 6.17% 

Me 1869.60 ±0.16 MeV 0.005% 

m c (m c ) 1.286 ± 0.053 GeV 4.10% 

m s (m c ) 0.122±g;g|g GeV 0.08% 

m b 4.651 ±0.054 GeV 0.03% 

V a , 0.2254 ±0.0013 0.03% 



Table 1: Hadronic (top), scale (mid) and experimental (bottom) uncertainties affect- 
ing the Standard Model calculation D + - D° lifetime difference in the MS scheme. 



yet, it is not an SU(3) breaking quantity and does not suffer from the large GIM 
suppression inherent to the latter. If the heavy-quark expansion is to be applicable, 
it will have to reproduce charm hadron lifetimes within the first few orders in l/m c , 
with possible new physics contributions expected to be small. For this, we think that 
meson lifetimes might be a good test for the inclusive framework in charm. 

Taking the lifetime ratio r (D + ) jr (D°) at leading order in the l/m c expansion, 
and comparing to the naive ratio of experimental averages [15], we find very good 
agreement: 



t(D° 



- = 2.536 ±0.019 

exp 



T (-^ + ) — 2 8 ± 1 5 ( hadronic ) +0 ' 3 ( scale ) _(_ q 2 ( 6X p) (7) 



MS 



t- (hadronic) +0.4 (scale) _i_ ri 9 (exp) 
^)pole -°- 9 ' 

Our Standard Model estimate is based on results from B physics [16] and includes 
QCD at next-to-leading order, renormalized in both MS and pole scheme. Uncertain- 
ties due to hadronic and experimental input and the variation of renormalization scale 
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Figure 2: The Standard Model estimate for r (D + ) /r (D°) at operator dimension six, 
as a function of the AC = 1 renormalization scale fii, in the pole (left) and MS scheme 
(right). The plots show the full result at NLO (solid) and the LO only (dashed). The 
error bands include all sources of uncertainty listed in Tab. [TJ we display them for a 
variation of the hadronic parameters 6%^ within ±0.1 (hatched) and ±0.05 (shaded). 
The thin horizontal line marks the ratio of the PDG lifetime averages. 



have been disclosed separately. For the MS result, we list their individual sources in 
Tab. (U to obtain the limits quoted in (J7J) above, the contributions shown in the Table 
have been added in quadrature. The dependence on the renormalization scale of the 
AC = 1 effective theory is illustrated in Fig. |2J Perturbation theory apparently 
becomes unreliable below about 1 GeV, yet it still seems to be under control at the 
charm threshold. Adding the QCD corrections significantly reduces the scale depen- 
dence. To quantify the scale uncertainty in (jTj), we varied the renormalization scale 
between 1 GeV and 2m c . The overall error is largely driven by hadronic uncertainties, 
entering through matrix elements of the dimension-6, AC = operators 

Q q = (cq) v _ A (qc) y _ A , 

Ql = (cg)s_ P (gc) s+p , 

r* = (cT a g) v „ A (gT a c) v „ A , ( } 

T| = (cT a q) s _ p (qT a c) s+P , 

where T a is the generator of colour SU(3). The meson state matrix elements of 
these operators enter r (D + ) jr (D°) in isospin-breaking combinations, conventionally 
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parametrized as [TU [18] 



<D+| 


Q u - 


-Q d | 


D+) 


<D+| 


Qs- 




D+) 


(D* 


| T u 


rpd 


D + ) 


(D+ 


1 T 1 ^ 

1 J s " 


rpd 1 

" J s 1 


D+) 


symmetry, (D° \Q U ' 


d \D° 


) = 



f 2 D M 2 D B 1 , 
f D M 2 D B 2 , 

9 9 ( 9 ) 

flM 2 D e x . 

(D + \ Q d,u \D + ). In vacuum saturation 
approximation (VSA), B\ — 1, B 2 — 1 + C(a s , Aqcd/wc), while e x and e 2 vanish. 
Non-factorizable corrections to the VSA values are of order 1/N C . At the b quark 
threshold, the hadronic parameters B\, B 2 , e±, e 2 are known from lattice QCD [T9] . 
Older studies rely on lattice HQET [20] and HQET sum rules [2T| . and disagree ma- 
terially with the former. All available non-perturbative studies, though, in agreement 
with the results for r (B + ) jr (B®) [16J, consistently indicate that |ei| and \e 2 \ are sig- 
nificantly smaller than 1/N C . The results from lattice QCD are, with the information 
made public, not easily convertible to charm. In our calculation we used the VSA 
values for all hadronic parameters and assigned errors 

(B u B 2 , ei , e 2 ) = (l ± y Nc , 1 ± y Nc , ± l/io, ± Vio) • (10) 

A lattice update for these quantities is urgently needed and will shrink the error on 
r (£>+) jr (D°) drastically. 

Our result (j7]) certainly supports some confidence in the OPE approach to charm 
physics. Both the closeness to the observed number and the reasonable behaviour 
of the QCD perturbation series do not show indications for a breakdown of heavy 
quark methods. Note, however, that this good agreement still might be by accident; 
barring significant uncertainties, we so far included only the leading order in l/m c , 
contributing through operators of dimension six. We expect sizeable corrections from 
the sub leading dimension seven. An estimate for these effects will be crucial to judge 
the convergence of the heavy-mass expansion. 



3 Charm mixing and SU(3) breaking 

Straightforward application of the heavy-quark expansion to the AC = 2 Hamilto- 
nian fails to predict mass and width splittings by orders of magnitude. At leading 
order in l/m c , charm mixing is predicted to be very slow due to severe GIM sup- 
pression. Interference among states in the same SU(3) multiplet almost cancels the 
D-D transition amplitude. The D-D width difference, for example, receives contri- 
butions from on-shell internal ss-, sd- and <i<i-pairs, see Fig. [3j In the limiting case 
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Figure 3: The l/m c - leading effect generating a D-D transition, contributing to the 
heavy-quark expansion through operators of dimension six. 

of exact SU(3), they cancel almost exactly (up to terms of order |Vrf,V^,| 2 ) as soon 
as the unitarity of the CKM matrix comes to work. Breaking of SU(3) enters the 
amplitude through a non-zero strange mass. To overcome the cancellation, one mass 
insertion per internal line is needed to break SU(3), and a second one to compensate 
the chirality flip. In the end, a factor of m^/m^. (in the CKM-leading part), remi- 
niscent of broken SU(3) symmetry, is suppressing the mixing rate. The interference 
will be prevented as soon as SU(3) breaking is introduced from sources other than 
the quark masses. If one of the internal momenta is less than Aqcd, the intermediate 
state couples to the meson's soft QCD substructure and feels the breaking of SU(3) 
in the hadron state. In the OPE picture, the effect is generated by operators of di- 
mension nine (Fig. HJ). Albeit suppressed by three additional powers of l/m c with 
respect to the leading dimension six, it is of lower order in m s and has been suspected 
to dominate the heavy-quark expansion [22]. Breaking of SU(3) symmetry in non- 
perturbative hadronic QCD dynamics can access the AC = 2 Hamiltonian through 
matrix elements of 6-quark operators containing a s- or <i-quark pair. To estimate 
the hadronic matrix elements, we assume that the fields from the intermediate state 
factorize from the operator structure. With the coupling to the sea quark background 
dominated by low energies, we model the hadronic s and d quark content with the 
vacuum condensate and evaluate the matrix elements of the factorized 2-quark struc- 
ture assuming vacuum saturation [23]. Corrections to factorization are of order 1/N C . 
The condensate of quark-antiquark pairs in the QCD vacuum is 

— -* *— = (Q\q(x)®q(0) IQ) = -J^ xl c (l-^£^ + ...Y (11) 
(qq) AN C \ A J 

for our purpose with q = d, s. The second term on the right-hand side corresponds 
to the first order term of a Taylor expansion in a small momentum flow through the 
hadron state. Higher orders in the expansion are suppressed by even more powers of 
m q and can be neglected. For the diagrams of Fig. H]to overcome chirality suppression, 
accordingly, we have to supply one mass insertion from the right-hand side of ( Hip , or 
from the perturbative quark propagation. Taken together, this amounts to a relief in 
SU(3) interference by one factor of m s /m c . If this relief is strong enough to outweigh 
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Figure 4: AC = 2 transitions coupling to the hadronic sea quark background in the 
intermediate state. Shown are the contributions to the D D width difference. One 
gluon is necessary to generate an on-shell intermediate state. 

the suppression by additional powers of l/m c (which isn't actually a so small number), 
diagrams of this topology would actually dominate the heavy-quark expansion. 

In our calculation of the contributions to Ar as shown in Fig. HI we can confirm 
the expectation from this power-counting. We found that, as expected, the correction 
to each individual SU(3) amplitude is small (in the percent range), as it should be if 
the l/m c expansion is to converge. Due to less pronounced residual flavour symmetry, 
however, it survives the cancellations affecting the l/m c - leading contribution when 
differently-flavoured intermediate states are added, and exceeds the latter by a factor 
of order ten. Likewise, the Standard Model prediction for the D D width difference 
is enhanced by a factor of order ten. A similar calculation has been done for the D D 
mass difference. Our results for x and y are 

x = (6 ±2) • 1(T 5 , 

(12) 

y = (8 ± 9) • 1(T 6 . 

The mixing rates still miss the experimental numbers ([I]) by a factor around 100 
and 1000, respectively. Nevertheless, the calculation shows that the SU(3) breaking 
hadronic state interactions in dimension nine indeed are the (so far) dominant effect 
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generating the D D decay width difference. The calculation also allows to extract 
the weak phase in mixing. Up to operator dimension nine, we find <fi = 1.81q ^. To 
estimate the physical amount of CP violation mixing might account for, assume that 
there is some mechanism, able to break also the yet remaining soft SU(3) interference, 
which originates from the continuing quark lines in Fig. HJ and from the mass- 
dependence of the QCD vacuum condensate ffTTj) . The effect will mainly work in 
amplitudes with real CKM-couplings and enhance the mixing rates while reducing 
the weak phase. If we further assume that it is powerful enough to push x and y to 
their average experimental values, we find that a weak phase in the range of one per 
mille up to (at most) one per cent might be accommodated. 

4 Summary and outlook 

Among the four neutral meson systems K, B, B s and D, charm is the most challenging 
to approach theoretically. In these proceedings, we gave a review on Standard Model 
calculations of the D D mixing rate. Weak interactions of mesons can be described 
in the framework of the heavy-quark expansion, which had great success predicting 
lifetimes, mixing and CP violation of B° and B^ 1 mesons with high accuracy. At 
Moriond 2012, LHCb presented the first measurement (> 5a) of the decay width 
difference Ar of neutral B s mesons, which had been expected to be most sensitive to 
violations of duality. The agreement with the Standard Model prediction is intriguing, 
a fortiori in that the energy released in B s — > D S D S , which dominantly generates the 
width splitting, is only around 1.4 GeV. Applied to charm, it is well known that 
the heavy-quark approach meets serious problems. Predicted D meson mixing rates 
typically fall short of the experimental averages by several orders of magnitude, which 
is commonly attributed to a failure of quark-level perturbation theory at the charm 
threshold. Some observations, though, argue that it might be possible to understand 
also the D-D width splitting in the Standard Model. Energy releases in D° decays 
contributing to Ar do not significantly deviate from the 1.4 GeV in B s . A closer 
inspection of the heavy-quark expansion showed that both NLO and l/m c corrections 
are reasonably small to allow for a convergence of the series. 

We presented first results on the D + - D° lifetime difference in the Standard Model, 
which are in good agreement with experiment. Corrections from QCD at text-to- 
leading-order are moderate and significantly reduce the scale dependence. Our result 
is affected by large hadronic uncertainties associated with the matrix elements of 
four AC = operators. Future lattice calculations could reduce these uncertainties 
drastically. We think that our result supports some confidence in the short- distance 
approach to charm. It will be essential, however, to quantify the effect of l/m c 
corrections which we expect to be sizeable. 

Neutral meson mixing is an SU(3) breaking observable; in the Standard Model, 
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the quark masses are the spurions of SU(3) symmetry breaking and mass and width 
differences scale with powers of m s /m c . Mixing is therefore is predicted to be very 
slow at leading order in the heavy-quark expansion. It had been expected that SU(3) 
breaking interactions with the soft hadronic background, contributing in higher orders 
of the l/m c -expansion, might be the dominant contributions to neutral charm mass 
and width splittings and explain the measured D D decay width difference without 
spoiling the overall convergence of the expansion. We investigated the contribution 
to Ar from 6-quark operators generating the meson-antimeson transition, which is 
the l/m c - leading contribution where the effect is at work. We find that it is able to 
enhance the width difference by a factor of oder 10. On the basis of our results, we 
argue that we can not exclude CP violation in mixing of up to one per cent. Future 
studies should address the contributions from 8-quark operators, where breaking of 
SU(3) might be even more significant. 
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